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ABSTRACT 

An analytical study was undertaken to determine the 
two-dimensional wall conduction effects in an internally 
finned, rotating heat pipe. The finite element method was 
employed to generate computer results for a copper condenser 
with triangular fins. Heat transfer rates were shown to be 
approximately seventy- five percent greater than that pre- 
dicted by an earlier, one-dimensional analysis. Heat transfer 
rates were found to be insensitive to rotational speed and 
fin half-angle. Due to numerical difficulties within the 
finite element program, no data was obtained for the finned, 
stainless steel condenser. 
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Constant multiplier of second degree term of parabolic 
temperature distribution equation^ °F/ft2 
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h External environment convective coefficient, 
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Other Symbols 
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Subscripts 
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I . INTRODUCTION 



A. THE ROTATING, NON -CAP ILLARY HEAT PIPE 

The rotating, non-capillary heat pipe is a self-contained 
device designed to transfer large quantities of heat in a 
rotating system. The heat pipe contains an evaporator, a 
condenser, and a working fluid as shown in Figure 1. The 
working fluid is vaporized by heat addition at the evaporator. 
The vapor flows toward the condenser end as a result of a 
vapor pressure differential. There, the vapor condenses on 
the inner wall, releasing its latent heat of vaporization. 

The condensate then returns to the evaporator where the cycle 
is repeated. 

The difference between the conventional heat pipe and the 
rotating, non-capillary heat pipe is the condensate return 
mechanism. A conventional heat pipe employs the capillary 
pumping action of an internal porous wick to transport the 
condensate to the evaporator. The disadvantage of this scheme 
is the resistance to condensate flow through the wick. Addi- 
tionally, the presence of nucleate boiling in the evaporator 
could create bubbles inside the wick structure, reducing the 
effectiveness of the pumping action [1]^. In the rotating 
heat pipe, rotation about its longitudinal axis establishes 
a centrifugal force field which, when combined with a slight 

^ Numbers in brackets indicate references listed in the 
Bibliography . 
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taper of the condenser section, forces the condensate toward 
the evaporator. At high rotation rates, the heat transfer 
rate is significantly increased over the rate achieved in 
conventional heat pipes of similar dimensions. 

B . BACKGROUND 

Previous analyses [2-9] of the rotating heat pipe have 
been conducted to determine its performance characteristics. 
These studies determined, among other things, that under 
normal circumstances, the heat transfer rate is limited pri- 
marily by the amount of heat which can be transferred through 
the condenser wall. 

In an effort to improve the condenser heat transfer per- 
formance, Schafer [8] estimated and compared the thermal 
resistances across the condensate film, the wall, and the 
outside convective film using Woodard's earlier data [9]. 

He showed that with the use of outside spray cooling and a 
thin copper wall, the controlling resistance is the inside 
film resistance. Schafer suggested that this resistance can 
be reduced (1) by the proper selection of working fluid, and 
(2) through the use of internal axial fins as sketched in 
Figure 2a. 

Machining axial fins into the inner surface of the con- 
denser would improve heat transfer performance through two 
mechanisms: a decrease in condensate film thickness near 

the apex of each fin, and an increase in inner wall surface 
area. However, with the addition of fins, the wall thermal 
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resistance would also increase due to the increase in 
effective wall thickness. Of course, the use of copper as 
the wall material, because of its high thermal conductivity, 
would minimize the added wall resistance. 

Schafer then developed an analytical model of the intern- 
ally finned condenser for the case of fins of triangular 
profile (Figures 2a, 2b). In the analysis, he assumed one- 
dimensional heat conduction through the wall and assumed a 
linearly decreasing temperature profile from the fin apex 
(point A) to the midpoint of the trough (point C) . Results 
from the analysis indicated significantly improved heat 
transfer rates for the finned condenser compared to the smooth 
condenser. He pointed out, however, that two-dimensional 
conduction effects may be significant so that a more detailed 
analysis of the heat conduction problem within the wall was 
necessary . 

C. OBJECTIVES 

The objectives of this study were, therefore, to (1) 
develop an analytical solution to account for two-dimensional 
conduction effects within the wall for the triangular fin 
profile, and (2) compare this solution with the solution 
derived by Schafer. 
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II. THEORETICAL PROGRAM 



A. REVIEW OF SCHAFER ANALYSIS [8] 

Schafer's study considered an internal fin geometry 
represented by Figure 2a. Since the sum of wall thickness 
and fin height was assumed to be much smaller than condenser 
radius, an average radius from any point on the axis to the 
condenser wall was taken to be the distance from the axis to 
a position midway between the base of the fin and the apex 
of the fin, r^ . The centrifugal acceleration term was then 
evaluated at this average radius, r^w^. 

The thin wall assumption also allowed the curvature of 
the wall to be neglected. The problem could then be analyzed 
using a planar conduction heat transfer model as shown in 
Figure 2b. Furthermore, because of repetitive geometry the 
problem could be reduced to the simpler geometry shown in 
Figure 3. Cartesian coordinates were then selected, as shown, 
to permit a simplified fin condensation analysis. 

Assumptions in the analysis were based upon the earlier 
work of Ballback [2] . These assumptions were similar to those 
used in the Nusselt analysis of film condensation on a wall. 
The more important of these assumptions are; 

(1) steady state operation, 

(2) film condensation, as opposed to dropwise condensation, 

(3) laminar condensate flow along both the fin and the 
trough, 

(4) static balance of forces within the condensate. 
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(5) one- dimensional conduction heat transfer through the 
film thickness (no convective heat transfer) , 

(6) no liquid-vapor interfacial shear stresses, and 

(7) no condensate subcooling. 

Additional assumptions in the analysis are noted later. 

The complete development of Schafer's finned condenser 
heat transfer theory is presented in his thesis [8] . However, 
since the fluid flow and heat transfer theory was applicable 
to the two-dimensional conduction heat transfer analysis, a 
summary of his work is presented here, prior to the develop- 
ment of the two-dimensional conduction model. 

1 . Fin Condensate Mass Flow 

A Nusselt analysis was performed for film condensation 
on the wall of the fin. In this analysis, the fluid momentum 
equation was constructed to account for the effects of condenser 
cone half- angle, (p , fin half -angle, a, and the centrifugal 
force field. The pressure differential term was neglected. 
Condensate flow on the fin was assumed to be in the z-direction 
only. 

An expression for fluid velocity, w(y,z), was then 
derived by substituting the shear stress equation for laminar 
flow, ^ , into the momentum equation and solving 

for w(y,z). Average velocity, w(z), was then found by inte- 
grating w(y,z) over y from y = 0 to y=5(z), the local film 
thickness. This average velocity was then placed in the fin 
mass flow equation to yield an expression for fluid mass flow 
rate along the fin per unit condenser length as a function 
of 6 (z) , 



14 



I 









» 



( 1 ) 



Mj(z) = 



P£ co^CR^ + X sin (|)) cos (i> cos a <S^(z) 



3y. 



2 . Trough Condensate Mass Flow Equation 

Fluid flow within the trough was assumed to be one- 
dimensional in the x-direction and of constant velocity across 
the width of the trough. Therefore, condensate flowing from 
the fins into the trough was assumed to accelerate immediately 
to the velocity of the trough fluid. Cross-sectional area of 
the trough condensate region was twice the area shown in 
Figure 4. 

The x-direction momentum equation for an infinitesimal 
fluid element in the trough was constructed with the pressure 
differential term neglected. The trough mass flow expression 
was then developed in a manner similar to the development of 
the fin mass flow equation (1) . The average trough velocity 
was derived and substituted into the trough mass flow equation 
to yield total mass flow through the trough region, 






p^w^(R +x sin sin <(> 






(e6*Cx) + sin a cosar5*(x)). 



C2) 



where the expression in the second parenthesis is the cross- 
sectional area of the trough condensate region. 

Total condensate mass flow was identical to trough 
condensate mass flow since all condensate eventually flowed 
to the trough. Therefore, once the total mass flow rate at 
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any position x was known, trough condensate thickness, 6*(x), 
could be determined using equation (2). 

3 . Condenser Heat Transfer Determination 

Schafer assumed zero heat flux boundary conditions on 
both sides of the wall section shown in Figure 4. Temperature 
at the fin apex was taken to be fluid saturation temperature 
since film thickness at this position was assumed to be zero. 

He then imposed a linearly decreasing temperature on 
the inner wall surface from the apex of the fin to the mid- 
channel position. The surface temperature expression became 

Ts - T^(z) = 32, C3) 

where 3 was an unknown constant. This choice of temperature 
distribution was arbitrary but allowed a simplified one- 
dimensional conduction analysis of the condenser wall. He 
then imposed the restriction that no heat transfer would 
occur through section BE of the wall (Figure 3) . All heat 
transfer through the fin surface, thus, flowed out of the 
wall through section FE. 

Heat transfer through the fin was equated with one- 
dimensional heat transfer through a rectangular geometry of 
the same width as section FE, and of height, t*, equal to t 
plus one-half the fin height (Figure 4) . Temperature of the 
upper surface of the rectangle was set equal to the tempera- 
ture at a point midway between positions A and B. The 
expression for heat transfer rate per unit condenser length 
through the wall and outer convective thermal resistance 
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of the fin became 



t CTj - ) - T„ ] 



• * 

Q^(z3 



t* 

^1“ 



‘out 



(4) 



It was then desired to develop an expression for 

• * 

Qj (z) in terms of the heat transfer rate through the fin 
condensate film. First, the equation for fin condensate 
film thickness was derived by differentiating (1) with respect 
to z and substituting the result into the differentiated 
equation for mass flow rate per unit condenser length. 



The expression produced was algebraically rearranged to form 
a differential equation in 6(z), 



6^(z)d6(z) = 



kf (Ts - Pf dz 



P£ 0)^ (R^ + X sin (J) ) h£g cos ({) cos a 



(6) 



Equation (3) was then substituted in (6) and the result 
integrated over z from z = 0 to z to give an equation for 
6(z), 



<Slz) = r- 



2 k£ 6 P£ 



p/ (R^ + X sin 4>)h £g cos 4) cos a 



1 V 

z (7) 



The value for fin heat transfer rate per unit condenser 
length was then evaluated using 

Q^(z*) = h£g M^(z*) (8) 
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where z* is defined in Figure 4. This expression was set 
equal to (4) to yield a transcendental equation in 3, 
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^O 

T* 



out 



1/ V4 3/ 

(2K^g)V‘* (p^oj) [ (R +x.sin <))) h^ cos (j) cos a] z* ^ 









(9) 



where was the width of the base of the fin. The expressions 
on the left and right hand sides of this equation were equiva- 
lent only when all boundary conditions were satisfied. The 
task was then to solve iteratively for 3 until equation (9) 
was satisfied. (z*) was then determined in equation (4) 
using the final iterated value of 3. 

Heat transfer rate from the trough region was evaluated 



using 



Q^Cx) 



(T - T ) 
^ s avg-^ 

6*(x) 




( 10 ) 



where calculated as an average temperature of the 

trough wall and e was the total trough width (Figure 4) . 

A flowchart description of the computer program used 
in the Schafer analysis is given in Figure 5. The axial 
length of the condenser was first divided into increments 
of equal length, Ax. Solution of the heat transfer problem 
was then started by assuming an initial trough condensate 
thickness, 6*^ = 0.00006752 feet. This initial value was 
based upon a mass flow analysis by Sparrow and Gregg of 
condensation on a rotating disk [10] . It was assumed in 
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Schafer's analysis that the effect on the Sparrow and Gregg 
solution of enclosing the disk within the condenser would 
be minimal. A study of the effects of a ten- fold increase 
or decrease from the above value of 6*^ indicated that the 
initial 6*^ assumption would not affect the heat transfer 
solution by more than ten percent. Therefore, the deter- 
mination of an exact value of initial trough condensate 
thickness was not crucial to the problem. 

Calculation of condensate mass flow contribution from 

the endplate then followed. Then, for each increment Ax, an 

3 z 

initial value of 3 was calculated by removing the term 

from the left side of equation (9) and solving directly for 
3. The full equation was then employed to solve iteratively 
for 3. Calculation of incremental fin heat transfer rate 
then followed by multiplying the value obtained in (4) by 
the incremental length. Ax. Incremental trough heat transfer 
rate was next evaluated by multiplying equation (lO) by Ax. 
Incremental trough and fin heat transfer rates were then 
summed to yield the total incremental heat transfer rate for 
the geometry of Figure 4. 

Total incremental mass flow rate was evaluated by 
dividing twice the sum of the incremental trough and fin heat 
transfer rates by the latent heat of vaporization of the fluid. 
As the incrementation progressed, incremental mass flow rates 
were added to mass flow from preceding increments to yield 
C^) • Equation (2) was then solved for 6*(x) with the 
aid of a polynomial rootfinder subroutine. The new value of 
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6*(x) was then used to begin the solution for the next 
increment. When the total length of the condenser had been 
traversed, total heat transfer rate was evaluated. 



^tot 



= 2 N 



fin 



100 

I 

i=l 



Ax[Q (z*). + 



Q,(x)i] 



B. TWO-DIMENSIONAL HEAT TRANSFER SOLUTION 



( 11 ) 



1 . Introduction 

In this study, two-dimensional conduction through the 
condenser wall was considered in the heat transfer analysis 
using a finite element solution. The finite element method 
was chosen over the method of finite differences because of 
considerations involving wall geometry and computer utiliza- 
tion time. Myers [11] has shown that when these methods are 
compared for the case of two-dimensional conduction in a 
rectangular block with heat generation, the finite element 
method more closely approximates the exact solution when the 
number of nodal points considered in each method are equal. 
This is particularly true when the nodal point meshes are 
coarse . 



To achieve the same accuracy as with the finite element 
method, the finite difference method must use a finer nodal 
mesh, especially for non-rectangular geometries. Additionally, 
when fins of different shape are considered, the mesh of the 
finite difference model must be completely reconstructed for 
each new shape. In the finite element program, however, a 
simple change in nodal coordinates is all that is necessary 
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to accommodate the change in fin shape. Since future studies 
may analyze the use of varied fin shapes with the rotating 
heat pipe, it was desired to incorporate the more flexible 
finite element method in the computer program to avoid having 
to rearrange the program for each fin shape. 

While the accuracy of the finite element method is 
superior, the finite difference method requires less computer 
utilization time for the same number of nodal points. If a 
finite difference model of finer mesh was chosen, the decrease 
in computer time due to the method might be offset by an in- 
crease in time required for the increased number of nodes. 

The overall improvement in computer utilization time might, 
therefore, be minimal. As a result, in this study, the choice 
of method was not based upon computer time requirements. Also, 
a finite element computer program [12] was available which 
could be adapted for the purposes of this analysis. For the 
reasons outlined above, therefore, the finite element method 
was selected. 

The steady state differential equation for conduction 
heat transfer is given in Arpaci's text [13] as 



3 , , 3T s ^ 3 . , 3T . ^ 

^ ( k ) + ^fk -rrrJ O* 

3y y 3y 3z ^ t a? 



"z 3 z 



( 12 ) 



Boundary conditions are normally of two kinds; those of 
specified temperature on a boundary. 



T = T^ , (boundary b only) 



( 13 ) 



and those of specified heat flux or convective condition 
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on the boundary, 




+ [h(T - T^) ] cijj = 0 . (boundary b only) 



(14) 



where H and i are direction cosines o£ the outward 
y z 

normal to the boundary's surface. 



2. The Finite Element Method 



Many references including [11], [14], and [15] develop 



the concept of using material elements of finite dimensions 
within a structure to model the differential equation and its 
boundary conditions. The steady-state field equation, (12) , 
may be discretized by considering that the temperature at any 
point within an element is a function of the temperatures of 
specified nodes on that element. 



field equation and boundary equations for each element may 
be mathematically rearranged and integrated to yield a matrix 
system of equations for that element. 



T = <N> { T }®, 



(15) 



where <(N^ is the row matrix shape function of position and 
{ T is the nodal temperature vector for the element. The 



[H]® { T = -{ , 



(16) 



where 




and 
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* (/^JN)h<N>dc)^ - MT„)^c^. 

Equation (18) is applied only over relevant boundaries. The 
element matrix equations are then assembled into a master 
matrix equation for the entire structure. Solution by 
numerical methods then follows. 

The convective boundary conditions were applied in 
the finite element solution of the finned condenser using a 
lumped convective condition formulation. Average nodal con- 
vective coefficients, h , were determined for each node 

avg 

subject to a convective environment. The nodal h values 

avg 

were then considered to be constant over the boundary region 
associated with the node. The method of determination of 
these coefficients is described in a later section. 

Equation (18) was, thus, reduced to the following 

form , 

{F} ® = ({h^,^„ T} T C19) 

b avg b avg <» b b 

This lumping procedure reduced the complexity of the convective 
condition calculations and did not seriously affect the values 
within the master matrix, [H] . 

Reference [15] presents a development of the 
isoparametric finite element concept. This concept allows 
the solution of equations (12) , (13) , and (14) for elements 
of non-uniform shape. The application of the method of finite 
elements is, therefore, extended to a wider range of physical 
problems . 
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The eight-nodal finite element mesh used in the 
two-dimensional, analysis of the axially finned condenser is 
shown in Figure 6. A computer program for the case of two- 
dimensional conduction was developed using a modified version 
of the three-dimensional, isoparametric finite element program 
assembled by Lew [12] . 

3 . Test Cases 

Before the finite element solution was applied to the 
condenser problem, a series of test cases were analyzed. 

These cases allowed a comparison of analytical solutions to 
results obtained using the finite element program. The two- 
dimensional heat conduction problems studied were 

(1) a rectangular block with specified surface temperature, 

(2) a rectangular block in a uniformly convective environ- 
ment, 

(3) a rectangular block with a triangular fin in a uniformly 
convective environment, 

(4) a rectangular block with a convex parabolic fin in a 
uniformly convective environment, and 

(5) a rectangular block with a concave parabolic fin in a 
uniformly convective environment. 

Results of the test cases are presented in Figures 7-11 with 
values obtained using the finite element program placed in 
brackets. For the first two cases, the heat transfer rates 
calculated with the aid of the finite element computer program 
agreed well with the exact, two-dimensional solutions. 

Although finite element heat transfer solutions were rela- 
tively accurate, nodal temperature values were slightly 
different from the exact values. This temperature discrepancy 
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led to numerical difficulties when the finite element solution 
was applied to the stainless steel condenser. These diffi- 
culties will be discussed more fully in Chapter III. 

For the finned block cases, analytical expressions 
were derived for heat transfer rates through extended surfaces 
using the one-dimensional fin analysis as described in [16] . 
Heat transfer rates calculated by the finite element program 
for this uniformly convective environment did not vary sig- 
nificantly from the one-dimensional results. Since test case 
results from the finite element program agreed favorably with 
the derived results, the program was assumed to be operating 
properly and was applied to the condenser heat transfer 
problem. 

4 . Condenser Boundary Condition Application 

The axially finned condenser problem was analyzed 
by considering the geometry of Figure 6. Boundary conditions 
for the problem were considered to be the same as those in 
the Schafer analysis with the exception of the inner surface 
temperature distribution. Instead of using a linearly decreas- 
ing temperature profile along the wall from the apex of the 
fin to the midpoint of the channel, a parabolic distribution 
was assumed and used to calculate average nodal convective 
coefficients . 

The Nusselt analysis of the fin condensation problem 
yielded a differential equation in 6(z), fin condensate 
thickness, which is identical to equation (6) used in Schafer's 
analysis . 
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Consistent with the eight-nodal finite element 
concept, a parabolic temperature distribution was then 
assumed to exist along the fin surface using the boundary 
condition, T = at z = 0, since 6(0) = 0. Therefore, 

T^(z) = a^z^ + bjZ + Tg , (20) 

or 



\ 

where a^ and b^ were unknown constants. 

When (21) was substituted into equation (6) and the 
result integrated from z = 0 to z , the expression for 6(z) 
became 



6(z) 




p£^ w^(R + X sin 4>) h£^ cos (}) cos a 



( 22 ) 



The local convective heat transfer coefficient, assuming 
one-dimensional heat conduction through the thin condensate 
film, was then 



h(z) 



fin 



^f 

^) 



= [■ 



k ^ 
^f 



w^(R + X sin <}) ) hr cos<|) cos a 



4p£ (- 



Z 3 



- b 









V4 

I 

'] 



(23) 



and was applicable from z = 0 to z = z* = z^ - 6*/cos a 
(Figure 6) . 



Across the width of the trough, from z = z^ to the 
mid-channel position. 



k£ 

trough 6*(x) 



(24) 
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6*(x) was considered constant across the trough width. Based 
upon a check of the wall temperatures generated from Schafer's 
computer program, the fluid temperature difference across the 
width of the trough was predicted to be small. Fluid thermal 
conductivity values were, therefore, assumed to remain 
relatively constant. As a result, was assumed to 

be constant across the trough width and was calculated using 
equation (24) and a representative average inner wall temper- 
ature to determine k£. 

An h(z) temperature distribution was then needed 

between z = z* and z = z . Since the thickness of condensate 

o 

film increased linearly from the value at z* to the value at 
z^, the convective coefficient was assumed to decrease linearly 
from h(z*)f.j^ to h(x)„ough- 

Other boundary conditions were then applied in 
accordance with the assumptions of Schafer, namely, zero heat 
flux through the sides of the geometry, T = at the fin 
apex, and constant convective coefficient, , along the 

outside surface of the wall. Boundary conditions for the 
problem were, thus, completely specified. 

In his expression for trough condensate cross-sectional 
area, Schafer assumed that the surfaces of the condensate 
flows from the fin and the trough were coincident at a position 
directly above the trough-fin wall juncture (Figure 4) . For 
small fin film thicknesses and other than small fin angles, 
this assumption was not accurate. In this study, it was 
assumed that fin condensate thickness was much smaller than 
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trough condensate thickness. The trough condensate region 
could, therefore, be more accurately described as shown in 
Figure 6. The modified trough mass flow equation ^became 

P£^ (R^ sin <}>)6*(x)^ sin <p 

" T (6*(x)e + 6*(x)2tana) . 

:>Uf 

(25) 

It should be noted that this equation differs from equation 
(2) only in the expression for cross-sectional area of the 
trough condensate region. The change in cross-sectional areas 
involved the term, 1/cos^ a, which for large fin angles, could 
make a significant difference in cross-sectional areas of the 
two equations. 

5 . Condenser Heat Transfer Solution 

Within the computer program, the heat transfer rate 
for the entire condenser was determined by dividing its axial 
length into one hundred equal length increments of length, ^x, 
solving for the incremental heat transfer rate using an iter- 
ative procedure, and summing those incremental rates to yield 
a total heat transfer rate. A flowchart description of the 
computer program used in this analysis is given in Figure 12. 

At each increment, constants a^ and b^ were determined 

using a Lagrangian interpolation fit of the temperatures 

/• 

calculated from the previous iterate at nodes 1, 3, and 6. 

For the first increment, a^ and b^ were initially calculated 
from an initial guess of the temperatures at those nodes. 
Convective coefficients at several locations along the fin 
were determined using equation (23) and the linearly decreasing 
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h(z) assumption from z = z* to z = z^. Next, average nodal 
values, hj^^g, and the average h o£ the left hand 

boundary region associated with node 6, (Figure 13), were then 
calculated using the position dependent values of hCz)^^^ in 
a Simpson's Rule approximation. trough then calculated 

from equation (24) using the 6* value calculated from the 
previous increment. This value was used as the average 
convective coefficient for nodes 9 and 12 and the right hand 
boundary region of node 6. The average convective coefficients 
calculated for the left and right hand boundary regions of 
node 6 were then weighted by their respective region lengths. 
The two values were then summed to yield the average convective 
coefficient for node 6, hg^^g. The remaining boundary condi- 
tions were then applied and incremental heat transfer rate 
per unit condenser length, , and nodal temperature values 
were generated from the finite element subroutines. 

If the resultant heat transfer rate converged to 
within the criteria of 0.05 percent change from the value of 
the previous iterate. 







< 0.0005 , 



(26) 



the incremental problem was deemed solved and set equal to 
the final iterate, Qj . If, on the other hand, the convergence 
criteria was not met, the newly calculated nodal temperature 
values were used to establish a new set of inner wall 
convective coefficients. After the remaining boundary 
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conditions had been applied, the finite element solution was 
repeated until heat transfer rate converged to within the 
acceptable 0.05 percent tolerance. Because Qj was initial- 
ized at zero for each increment, the iteration procedure was 
repeated at least twice for every increment. 

Incremental mass flow rate was then calculated using 



A M 



tot. 




(27) 



This value was then added to the mass flow rate from previous 
increments. The subsequent interval's trough condensate 
thickness was then calculated using a polynomial root-finder 
subroutine to solve equation (25) . The incrementation 
continued until the entire length of the condenser had been 
traversed. Incremental heat transfer rates were than summed 
and a total heat transfer rate calculated from (3) . 

100 • 

•5tot “ 2 Nf (Q. AX) , (28) 

where was the number of fins along the inner wall 

circumference . 

To start the incrementation, two initial assumptions 

were made. The initial increment's value of 5* = 6* was 

1 

taken from the analysis of Schafer as 0.00006752 feet. For 
reasons explained in the previous sub-chapter, a large variance 
in the accuracy of this value would not significantly affect 
the heat transfer results obtained from the computer program. 
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Additionally, at the first increment, initial guesses 
of fin temperatures at nodes 1, 3, and 6 were made in order 
to begin the solution. As an example, the temperatures used 
to begin the solution for the copper condenser with fluid 
saturation temperature of 200®F were = 200®F, = 195°F, 

and T = 191°F. 

6 

For the copper condenser, the use of widely variant 
initial nodal temperature values, = 200°F, = 191°F, 

and T = 182®F had little effect on convergence to the heat 
transfer solution for the first increment. However, as will 
be discussed in Chapter III, the initial temperature guess, 
coupled with the slight inaccuracy of the finite element 
nodal temperature solution produced numerical difficulties 
within the computer program for the stainless steel condenser. 
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III. DISCUSSION OF RESULTS 



Computer results for the finned copper condenser were 
generated for three outside convective coefficients, three 
fin angles, and five rotation rates. The working fluid used 
throughout the analysis was water. Representative results 
are presented in Figures 14 and 16, 17, 18. 

Figure 14 compares the heat transfer results of this 
study with the results of the analysis of Schafer for given 
operating conditions. It can be seen that the two-dimensional 
conduction analysis allowing a parabolic temperature distri- 
bution along the fin predicted heat transfer rates for the 
finned condenser which are approximately seventy-five percent 
greater than the rates predicted with the one-dimensional 
model with linear temperature profile along the wall. This 
rather significant change was unexpected, and may be some- 
what in error due to numerical difficulties in the method. 
These numerical difficulties arose due to the boundary 
condition at node 1. The finite element boundary conditions 
are established for each node such that the condition is 
assumed to exist over the entire boundary region of the node. 
Therefore, the zero flux boundary condition could not be 
applied over the whole length of the left and right boundaries 
of the condenser geometry. This problem existed for nodes 1, 
4, 11, and 12, but was most critical for node 1 since it was 
estimated that twenty to thirty percent of the total heat 
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flux through the inner wall passed through the region 
associated with this node. 

An error of fifty percent in the calculation of node 1 
heat transfer rate might then give a total heat transfer 
error of as high as fifteen percent. Two possible solutions 
to the boundary condition problem merit attention. The first 
is to increase the number of elements used in the solution of 
the geometry of Figure 6. However, a boundary region associ- 
ated with node 1 would still exist on the left hand side of 
the geometry although reduced in size. The effects of the 
boundary condition problem would, therefore, be reduced but 
not eliminated. 

The second solution would be to consider a wall geometry 
as shown in Figure 15. The zero flux boundary condition 
would then be forced along the boundary in question due to 
the symmetry of the problem. The use of one of these alter- 
nate solution methods might contain the added benefit of 
eliminating the numerical difficulties associated with the 
stainless steel condenser solution. 

The heat transfer rate solution is generally affected in 
the same manner as with the analysis of Schafer when h 

out 

0 ), and a were varied. In Figure 16 the influence of outside 
convective coefficient, , is shown. Figure 17 illustrates 

how the rate of heat transfer was relatively unaffected as 
fin half angle, a, was varied. In addition, heat transfer 
rate was found not to vary significantly with rotation rate 
as shown in Figure 18. This result concurred with the results 
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obtained earlier by Schafer [8] . Schafer suggested that the 
invariance in heat transfer rate with the rotation rate was 
due to the cancelling effects of decreasing film thickness 
and decreasing temperature difference across the film 
thickness as rotation rate is increased. 

Because of numerical difficulties in the program, no 
results were generated for the stainless steel condenser. 

With the solution of the first and second test cases, the 
nodal temperatures calculated were accurate only to within 
two. or three degrees. In the condenser problem, this slight 
inaccuracy resulted in nodal temperatures along the fin wall 
which were slightly above the fluid saturation temperature. 
Computer output indicated that the wall temperature at node 3 
was approximately 0.1°F above the saturation temperature when 
Tg = 200°F. As a result, the quantity within the brackets of 
equation (23) was a negative value. Error messages were 
generated when the computer attempted to evaluate the fourth 
root of this negative number. 
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IV. CONCLUSIONS 



1. The addition of internal axial fins to the condenser 
surface of a rotating heat pipe significantly improves heat 
transfer performance, 

2. Analysis of the condenser heat transfer problem using 
a two-dimensional wall conduction model and a parabolic 
temperature distribution along the fin yielded heat transfer 
rates about seventy-five percent above those predicted by the 
one-dimensional, linear temperature distribution analysis of 
Schafer [8] . 

3. For the geometrical cases studied in this thesis, the 
heat transfer rates calculated from the two-dimensional 
analysis are insensitive to changes in rotational speed and 
fin half-angle. 

4. Boundary condition problems produced with the finite 
element solution for the copper condenser were the probable 
cause of solution inaccuracy on the order of fifteen percent. 

5. Slight inaccuracy from the finite element calculation 
of nodal temperature values led to numerical difficulties in 
the heat transfer solution of the stainless steel condenser. 
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V. RECOMMENDATIONS FOR FURTHER STUDY 



It is suggested that the following topics be considered 
for future study of the internally finned, rotating heat 
pipe : 

1. Re-examination of the finite element analysis using 
either an increased number of elements for the geometry 
considered in this study or the mirror image geometry pre- 
sented in Figure 15. 

2. Use of twelve-nodal elements in the finite element 
analysis to allow the assumption of a cubic inner wall temper- 
ature profile in order to more closely approximate the actual 
temperature distribution. 

3. Analysis of film condensation on a variety of fin 
shapes, and application of the solution to the two-dimensional 
model to determine the shape permitting optimum heat transfer 
rate. 

4. An experimental program in order to compare theoretical 
and experimental heat transfer rates of the internally finned 
condenser. 
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Figure 1. Rotating Heat Pipe Components 
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Figure 2a. Internally Finned Condenser Geometry 
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Figure 2b. Planar Conduction Simplification 



38 






] 

i 

1 




♦ 



y 




Figure 3. Condenser Geometry Considered in Analysis 
of Schafer [8] 
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<D = 1° 



Figure 4. 



Condenser Geometry Considered in Analysis 
of Schafer [8] 



40 



i 

i 

i 



I 




Figure 5. Computer Program Flowchart for the Analysis 
of Schafer [_8] 
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Figure 6, Condenser Geometry Considered in 
Two-Dimensional Analysis 
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Figure 7. Test Case (1) 
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Figu re 8. Test Case (2) 
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Figure 9. 



Test Case (3) 
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Figure 10. Test Case (4) 
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Figure 11. Test Case (5) 
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Figure 12. Computer Program Flowchart of 

Two-Dimensional Conduction Analysis 
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Figure 13. Node 6 Boundary Regions 
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Condenser Heat Transfer Rate, BTU/Hr (x 10'^) 




Saturation Temperature, °F 



Figure 14. Condenser Heat Transfer Rate vs. 
Saturation Temperature 
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Figure 15. 



Alternate Condenser Wall Geometry 
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Condenser Heat Transfer Rate, BTU/Hr (x 10"*) 




Figure 16. Condenser Heat Transfer Rate vs. 
Saturation Temperature 
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Condenser Heat Transfer Rate, BTU/Hr (x 10"^) 




Figure 17. 



Condenser Heat Transfer Rate vs. 
Saturation Temperature 



Condenser Heat Transfer Rate, BTU/Hr (x 10"^) 



h 1000 Hr-ft2- °F 
a = 30° 




Figure 18. Condenser Heat Transfer Rate vs. 
Rotation Rate 
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